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Abstract. We present constraints on cosmological and astrophysical parameters from high- 
resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cos- 
mology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of 
primary cosmological and secondary foreground parameters is fit to the map power spec- 
tra and lensing deflection power spectrum, including contributions from both the thermal 
Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson 
and correlated anisotropy from unresolved infrared sources, radio sources, and the correla- 
tion between the tSZ effect and infrared sources. The power 2-k of the thermal SZ power 
spectrum at 148 GHz is measured to be 3.4T1.4 /rK 2 at £ — 3000, while the corresponding am- 
plitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 /jK 2 . 
Combining ACT power spectra with the WMAP 7-year temperature and polarization power 
spectra, we find excellent consistency with the LCDM model. We constrain the number of 
effective relativistic degrees of freedom in the early universe to be iV e ff = 2.79±0.56, in agree- 
ment with the canonical value of lV e ff = 3.046 for three massless neutrinos. We constrain 
the sum of the neutrino masses to be E m u < 0.39 eV at 95% confidence when combining 
ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain 
the amount of primordial helium to be Y p — 0.225 ± 0.034, and measure no variation in the 
fine structure constant a since recombination, with a/ao — 1.004 ± 0.005. We also find no 
evidence for any running of the scalar spectral index, dn s /d\nk — — 0.004 ± 0.012. 

Keywords: CMBR experiments, Sunyaev-Zeldovich effect, physics of the early universe, 
cosmological parameters from CMBR 
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1 Introduction 


Studies of the cosmic microwave background (CMB) have dramatically progressed over the 
past two decades [e.g., 7, 19, 33, 35, 36, 45, 73, 75, 99, 105, 112, 115, 130, 136, 139, 154, 157, 
168, 169]. The current ACDM cosmological model provides an excellent fit to the CMB data 
across a wide range of angular scales, and is supported by complementary observations of 
large scale structure [e.g., 92, 152], Baryon Acoustic Oscillations [BAO e.g., 6, 24, 27], Type 
la supernovae [e.g., 38, 87, 102], galaxy cluster measurements [e.g., 120, 155, 180, 189] and 
observations of gravitational lensing [e.g., 63, 86, 103, 125, 156, 175]. While the CMB on 
angular scales of greater than 0.3° has been definitively measured by the Wilkinson Microwave 
Anisotropy Probe [WMAP, 17], a wealth of information in the CMB on smaller angular 
scales continues to be probed with ever-increasing precision [e.g., 44, 85, 101, 113, 145- 
147, 150, 151, 164]. As this paper was being finalized the SPT collaboration released a 
new set of papers [95, 173] and the WMAP team released its final 9-year results [17, 90]. 
Our analysis does not incorporate these results although in a few places we make direct 
comparisons. While this paper was under peer review, the Planck satellite released its first 
cosmological results. We leave a final combined ACT+Planck analysis to future work. 

The Atacama Cosmology Telescope (ACT) complements measurements from WMAP 
by observing from l ~ 300 to £ = 10000. This widens the range of data available to constrain 
both cosmological parameters through the Silk damping tail of the primary CMB [166] and 
the residual power from secondary sources between us and the surface of last scattering. These 
sources include galaxy clusters, which are detectable at microwave frequencies through the 
Sunyaev-Zel’dovich (SZ) effect [174, 195]. The thermal SZ (tSZ) effect describes the spectral 
distortion due to the inverse-Compton scattering of CMB photons to higher frequencies 
by the hot gas in clusters, while the kinematic SZ (kSZ) effect measures the corresponding 
temperature shift due to the bulk peculiar motion of the clusters. While these effects produce 
a diffuse signal at small scales from unresolved clusters, their influence has recently been 
detected directly through the cross-correlation of the ACT temperature maps with other 
tracers [76, 77]. In addition to emission from clusters via the SZ effect, radio galaxies and 
dusty star- forming galaxies also contribute to the power on small scales, and indeed dominate 
the cosmological signal for multipoles £ > 3000 and frequencies v > 90 GHz. Gravitational 
lensing by structures along the line of sight also generates a microwave background signal, 
and distorts the primordial CMB. 

This paper forms part of a set of papers presenting the 3-year analysis of the ACT 
data; the ACT temperature and deflection power spectra are presented in [42], while [54] 
presents the likelihood used in this analysis. [82] presents a catalogue of SZ-detected clusters 
from the ACT data, and interprets them. This paper contains the parameter estimation of 
both primary (cosmological) and secondary (foreground) parameters. We outline the data 
used in this analysis in section 2 and describe the methodology and likelihood in section 3. 
We present constraints on primary cosmological parameters in section 4 and on secondary 
parameters in section 5. We conclude in section 6, after which we provide an appendix of 
analysis tests. 

2 Data 

This paper presents results from a combination of observations at two frequencies, 148 GHz 
and 218 GHz, of multiple fields taken over three years. The southern fields (ACT-S) varied 
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over different seasons, with the 2008 season containing a 292 deg 2 patch and the 2009 and 
2010 seasons focusing on a smaller 146 deg 2 footprint. Equatorial data (ACT-E) were taken 
only over the 2009 and 2010 seasons; we use a 300 deg 2 patch for the 3-year analysis. We 
follow a similar procedure for going from maps to the temperature power spectra as was used 
in [44], using the power spectrum estimation procedure presented in [41]. 

The data and map-making procedure are described in [56] ; the power spectrum method 
and systematic tests are presented in [42] . Using the spectrum presented in [42] , we have con- 
structed two likelihoods; these likelihoods that are also presented in [54]. The multi-frequency 
likelihood parameterizes the foreground emission using additional parameters which we list in 
section 3, while the CMB-only likelihood marginalizes over these foregrounds. For the anal- 
ysis presented here we use the multi-frequency likelihood; however we show in appendix D 
that the two likelihoods give equivalent results. 

The data used in this analysis are the multi-frequency temperature spectra estimated 
from the ACT maps with combined spectra from both ACT-S and ACT-E including their 
covariance. The spectra are presented in [42] . For the results presented in figure 4 we use the 
marginalized CMB-only likelihood. In addition, we include the measurement of the power 
spectrum of the lensing deflection angle, £ 2 /4 Ci dd [42], 

The lensing spectrum is estimated from the ACT temperature maps using an optimal 
quadratic estimator [98]. Only data from the ACT equatorial patches are used to measure 
the deflection power, as the signal-to-noise of the southern patch was much lower than that 
of the equatorial data. The covariance between the lensing power spectrum and the tem- 
perature power spectrum is small. When adding in the deflection data, we will use the 
abbreviation ‘ACTDefT. 

The maps are cross-correlated with WMAP7 maps in order to obtain a calibration factor 
in multipole space. The cross-correlation calibration method is described in [72]; details of 
the calibration of the ACT 3-year data are given in [42], 

2.1 Beam and calibration errors 

Understanding the beam profiles is essential for interpreting the high-t 1 aspects of the power 
spectrum. At £ — 3000 and 148 GHz, the window function is ~ 0.5 its value at £ — 200. The 
beams are estimated independently for each array and season [56] from observations of Saturn 
and Uranus. The beams vary slightly with season due to changes in the telescope focus. We 
include a contribution in the likelihood to the full covariance matrix from the covariance of 
the beams for each season. The beam error includes contributions from the uncertainty in the 
pointing variation of the telescope. At £ — 700, the pivot point for the beam and calibration 
uncertainties, the effective calibration error is 2% for the 148 GHz maps and 2.6% for the 
218 GHz maps. Different seasons have somewhat different calibration uncertainties and so 
these numbers should be considered as representative of the effective combined calibration. 
While the absolute calibration is performed from cross-correlations to WMAP data, the 
telescope pointing solution, beams, and detector responsivity are characterized independently 
in each observing season and thus the calibration uncertainties of ACT-E and ACT-S are 
relatively independent. 

In the cosmological analysis, we apply a calibration prior for the 148 GHz spectra ob- 
tained from the WMAP- ACT cross-correlation calibration procedure described above. In the 
chains, we allow for a small error in the overall calibration of the spectrum by marginalizing 
over independent calibration factors for the south and equatorial spectra, at both 148 GHz 
and 218 GHz. This extra calibration allows for the overall spectra to adjust themselves at 
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the 1% level and has a negligible effect on the cosmological parameters, as discussed in ap- 
pendix C. Similarly, in the same appendix, we test for the dependence of the cosmological 
parameters on the beam’s assumed uncertainty, and find that beam error has a negligible 
effect on the parameters of interest. 

2.2 Additional data 

We use temperature and polarization data from the seven year data release of the WMAP 
satellite [109, 114] in addition to the measurements of the microwave temperature from 
ACT. We include measurements of the Baryon Acoustic Oscillations (BAO) from the Six- 
degree Field Galaxy Redshift Survey [6dFGRS, 20] and the Sloan Digital Sky Survey Data 
Releases 7 [SDSS DR7, 142] and 9 [SDSS DR9, 6], measured at redshifts: z — 0.106,0.35 
and 0.57. In addition, we supplement our data with a measurement of the Hubble constant 
of H 0 = 73.8 ± 2.4 km s -1 Mpc -1 [153], although [62] find H 0 = 74.3 ± 2.1 km s -1 Mpc -1 . 
In some cases we also include a prior on erg from skewness measurements of the tSZ effect 
from ACT [192], 

Unless explicitly specified, the ACT 3-year data are combined with WMAP7 data. 
For some model constraints, such as those on N e g and the secondary parameters, we also 
include the ‘low-U and ‘high-t 1 ’ spectrum measurements from SPT [101, 151], following the 
prescription in [54]. We show the ACT and WMAP7 data in figure 1. The best-fit model for 
the various frequency components is shown in figure 2. 


3 Methodology 

We use Markov Chain Monte Carlo (MCMC) methods to determine parameters associated 
with a variety of models. The basic cosmological ACDM model consists of 6 parameters 
describing a flat universe. These include the physical baryon density Hf,h 2 (where h — 
Hq/100 km s _1 Mpc^ 1 is the dimensionless Hubble parameter), cold dark matter (CDM) 
density fi c /i 2 , and 9a — r a /dA, the ratio of the acoustic horizon to the angular diameter 
distance at decoupling. This parameter is sensitive to the dark energy density, Ha, but less 
degenerate with other parameters [110]. The value of Ha is then a derived parameter. We 
assume the primordial perturbations to be scalar, adiabatic, and Gaussian and parametrize 
them via a spectral tilt n s , and amplitude A^, defined at pivot scale ko — 0.002 Mpc -1 . 
We assume that the universe transitioned from a neutral to an ionized state over a small 
redshift range, A 2 = 0.5, with optical depth r. The reionization history of the universe can 
be probed by small-scale CMB measurements through the impact of reionization on the kSZ 
effect [68, 106, 137, 194], although care must be taken to allow for correlations between the 
tSZ effect and the microwave emission from unresolved dusty galaxies [3, 129]. 

We express the ACDM set of parameters as 


{flb/i 2 , H c h 2 , 8a, A^j, n s , t } . (3.1) 

In this ACDM model, the number of effective relativistic degrees of freedom is assumed to 
be N e g — 3.046, with the abundance of primordial helium fixed at Y p — 0.24. The likelihood 
used in the ACT analysis is described in [54], which we briefly summarize here. We fit a 
model of secondary emission to the ACT multi-frequency power spectra that includes an 
additional nine parameters when considering ACT data in combination with WMAP7 data. 
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Figure 1. Data used in this cosmological analysis. The data from the WMAP 7-year data release [109, 
114] are combined with the ACT data. In this figure we show a weighted co-added spectrum from the 
equatorial and southern patches at 148 GHz. The full ACT likelihood, however, considers independent 
southern and equatorial spectra for both 148 GHz and 218 GHz [given in table 3 of 42, and online on 
LAMBDA] and their cross-frequency spectrum, which are shown in figure 2. The solid line indicates 
the best-fit cosmological model including foreground emission, while the dashed line shows the best-fit 
primordial CMB spectrum. This binning was selected for the common analysis of the equatorial and 
southern data. See [42] for alternative binnings. 


For the 148 GHz data we use modes 500 < t < 10000, while for the 218 GHz data we restrict 
ourselves to 1500 < £ < 10000. The theoretical spectrum for frequency bands i and j is 


°l 


_ k?CMB 

— 


+ B 


sec,ij 
i > 


(3.2) 


where B n — £{£ + l)Ci/2n and is the lensed primary CMB power spectrum. The 

secondary spectra components are modeled as 


R sec,ij 


Bf ' Z,ij + Bf [ Z ’ ij 

+ + Bf z ~ cm ^ + + H £ Gal ’ ij 


(3.3) 


with contributions from the tSZ and kSZ effects, CIB sources, the cross-correlation between 
the tSZ and CIB signals (tSZ-CIB), radio galaxies (rad), and residual Galactic dust (Gal). 
In addition to the six primary parameters, we add the following nine parameters 


{atsz, «ksz , ( Ip , & c , CLgsi dge: Pci &S1 £} • 


(3.4) 
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Symbol 

Description 

~ a 

«tsz 

Thermal Sunyaev-Zel’dovich (SZ) power amplitude. 

,, a 

Oksz 

Kinematic SZ power amplitude. 

a p a 

Poisson Cosmic Infrared Background (CIB) power amplitude. 

a c a 

Clustered CIB power amplitude. 

a gs a 

Residual galactic emission amplitude for the ACT-S spectrum. 

a ge a 

Residual galactic emission amplitude for the ACT-E spectrum. 

Pc 

Emissivity index of the clustered CIB power. 

a s a 

Radio Poisson power amplitude. 

C 

tSZ-CIB correlation amplitude. 


a Amplitudes relative to templates normalized to t(l + l)Ce/2ir = 1/rK 2 at = 3000 
and 150 GHz. See [54] for a complete description. 


Table 1. Secondary parameters in the ACT 3-year foreground model. 


The parameter atsz parameterizes the amplitude of the tSZ power; aksz the kSZ amplitude; 
a p and a c the Poisson and clustered Cosmic Infrared Background (CIB) power and a gs ,a ge 
model the residual Galactic dust anisotropy in the southern and equatorial survey regions. 
All the parameter amplitudes are dimensionless, and are defined for a template spectrum 
normalized to 1 /iK 2 at — 3000, and frequency v — 150 GHz. The frequency dependence 
of the correlated and Poisson CIB power is given in flux density units by the product of 
modified blackbodies with effective temperature 9.7K and emissivity index /3 C , following [4], 
and as described in Equations (8) and (9) of [54]. The radio source power has an amplitude 
a s and a spectral index a s fixed to —0.5, based on the assumption that the index obtained 
for brighter sources from ACT and SPT source catalogs [121, 188] holds for fainter sources. 

The clustered and Poisson (both CIB and radio) templates vary with scale as £(£ + 
l)Ce/2n tx £ 0,8 and oc £ 2 , respectively. We allow for a correlation between the tSZ effect 
and CIB sources, with scale dependence given by the template calculated by [3], and a 
frequency-independent correlation coefficient, £, which is defined in Equation (11) of [54], 
and is restricted to lie in the range 0 < £ < 0.2. The parameters in our foreground model 
are summarized in table 1. 

The spatially variable Galactic emission has been masked, leaving only a small residual 
component. To determine the amplitude and spectrum of the residual component we cross 
correlate with the IRIS map [131] as described in detail in [42], We then marginalize over this 
dust component in the likelihood with separate amplitudes in the equatorial and southern 
regions. There is roughly double the amount of dust in the equatorial region than in the 
south. In addition, there are some bright clouds. Investigations were performed to identify 
possible residual dusty clouds that were missed by our treatment, but no clear evidence for 
them was found. The likelihood marginalization described in [54] was found to be the most 
general and parsimonious treatment. 

We mask all sources above a detection threshold of 15 mJy. The source detection algo- 
rithm is described in [121]. The 148 and 218 GHz source samples in the southern map will be 
presented in [122], and source catalogs for the full data set will be presented in [67]. Sources 
are masked to a flux level of 6.4 mJy in the SPT analysis presented in [151], resulting in a 
source amplitude at 150 GHz of a' s — 1.3 ± 0.2. Sources are masked to a level of 50mJy 
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Figure 2. The ACT multi-frequency data. Each solid line indicates the best-fit total model for 
each cross spectrum, while the dashed lines show the radio power; the dot-dashed line indicates the 
power from CIB sources (both Poisson and correlated CIB); the gray lines show the total SZ power. 
The solid and dotted gray lines are for those spectra including 218 GHz and consist of only the kSZ 
contribution, while the dashed gray line is for the 148x148 GHz spectrum. The dotted lines show the 
(negative of) the (SZ-CIB) correlation. In all cases the models are for the best fit in a ACDM scenario. 


in [101]. The estimated difference in unresolved radio source power is 9.2 /rK 2 . Hence we 
include a separate amplitude for radio emission measured by the South Pole telescope, a' s , 
when including the SPT data in our analysis, whereas the other parameters in the secondary 
model, apart from the Galactic dust parameters, are common to both data sets. Thus we 
first subtract an amplitude of radio Poisson power of 9.2 fiK 2 from the [101] ‘low-^’ spectrum. 

For the parameter analysis, we use the publicly available CosmoMC code, which includes 
version 1.5 of the Recfast code [5, 37, 158, 159, 176, 193]. 

We combine measurements of the three independent lensed cross spectra: 148x148 GHz, 
148x218 GHz and 218x218 GHz made from the ACT-S and ACT-E fields described in 
section 2. In addition, we use the measurement of the lensing deflection field by ACT, 
presented in [42], We investigate four types of fits in this analysis: 

• We apply the full ACT likelihood to the data in combination with WMAP7 data and 
obtain constraints on both the primary and secondary parameters, with a full 15, 16 
or 17 parameter model (and four additional calibration nuisance parameters). 

• We estimate C i bandpowers marginalized over the secondary foregrounds, and obtain 
constraints on primary parameters based on these marginalized bandpowers. 

• We combine ACT and SPT to check for consistency between these small-scale 
experiments. 

• We consider the ACT data alone without WMAP7 data, while fixing (or placing priors 
on) the spectral index n s and the optical depth r. 
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Figure 3. Parameter constraints on the ACDM model for the combined ACT and WMAP 7-year 
power spectra. For each parameter, the full likelihood is marginalized over the other parameters. 
Dotted curves for the primary parameters are for WMAP7 data only. The solid curves are for 
ACT plus WMAP7, and the dashed lines add priors from measurements of the baryon acoustic 
peak [6, 20, 142] and a direct measurement of H 0 [153]. The constraints are summarized in tables 2, 4 
and 5. Since the £ parameter is unconstrained and set only by the prior, 0 < £ < 0.2, it is not shown 
in the tables. 


The models obtained from fitting only the ACT-S and ACT-E spectra are consistent 
with the models fit to the combined data, the best-fit spectra from each region agreeing to 
within 4%. Moreover, the best-fit theoretical spectrum from the current 3-year ACT data 
agrees with the spectrum derived from the 1-year ACT data to the 1% level. We discuss the 
consistency of the ACT spectrum in appendix A. 

4 Constraints on primordial parameters 

The ACDM model continues to fit the ACT data well, when combined with the independent 
WMAP7 data. Figure 3 illustrates the constraints on the ACDM model with the additional 
secondary parameters for the ACT+WMAP7 data combination. In addition, we plot the 
constraints from the WMAP7 power spectrum alone. ACT extends the angular range mea- 
sured by WMAP, but the parameters from the joint fit are consistent with those from WMAP 
alone. In addition, the plot shows that the six parameters are robust to the presence of low 
levels of foreground emission that can be identified and extracted by ACT because of its 
higher resolution and different frequency coverage. 

We start by constraining the parameters in the ACDM model. Our constraint on the 
scalar spectral index is n s — 0.971 ± 0.009, using ACT data in combination with WMAP7 
data, BAO and Hq measurements. We discuss the constraints on n s in section 4.5. In 
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addition, we improve the constraints on the baryon density using only CMB data, to 

ilfji 2 — 0.0225 ± 0.00047, which is due to the fact that the ACT spectrum places tight limits 
on the positions of the higher order peaks below t ~ 3000. The models are summarized in 
tables 2 and 4 for the ACT data in combination with WMAP7. Table 5 shows the constraints 
when adding BAO and Ho measurements to the ACT power spectrum data. In appendix D 
we show that the ACDM parameters derived from the CMB-only likelihood agree with the 
full likelihood to within 0.1a. 1 

The greatest power of ACT comes when quantifying models beyond the standard cosmo- 
logical model because the temperature power spectrum contains little additional information 
on the simple ACDM parameters at angular scales smaller than the third acoustic peak, 
i ~ 800 [e.g., 110]. The damping of the higher-order acoustic peaks relative to the base- 
line model [108] provides constraints on a variety of non-standard models. The ACT 1-year 
spectrum data [44] showed a slight excess of damping at small scales relative to the baseline 
model. Evidence for this slight excess is not seen in the ACT 3-year data set. More data 
resulted in a spectrum with smaller error bars, providing tighter constraints on parameters 
such as the baryon density, while the best-fit theoretical spectra are consistent between the 
two results at the 1% level. The consistency between the results presented here and those 
presented in [55] is discussed in section D of the appendix. 

Various authors have explored possible models which lead to excess damping of the Silk 
damping tail [e.g., 1, 30, 31, 58, 59, 64, 74, 81, 127, 128]. In this analysis we broaden our 
standard picture with other parameters and interpret the damping tail of the ACT data in 
this context. 

4.1 ACT data alone 

In figure 4, we compare the ACT and WMAP7 data both separately and together with the 
constraints on the ACDM model. This provides an important cross-check of the consistency 
of the two data sets. A key parameter which is primarily constrained with polarization data 
on the largest scales, such as those probed by WMAP, is the optical depth r. Hence, when 
considering the constraints from the ACT data alone, without including WMAP7, we impose 
a prior on the optical depth, r = 0.089 ± 0.015. In figure 4 we show two cases, one in which 
the scalar spectral index is fixed at n s — 0.967 and one in which n s is allowed to vary freely. 
Fixing n s tightens the bound on the amplitude of fluctuations, while other parameters are 
largely insensitive to the effect. The agreement shows that the same model that describes 
the WMAP7 data for f < 1000 independently fits the damping tail measurement from ACT 
of £ > 500. The same behavior is observed in [173]. 

4.2 Effective number of relativistic species 

The standard cosmological model has three neutrino species, all of which have negligible 
mass and contribute to N e g, the effective number of relativistic species at recombination. 2 
Precision electroweak measurements place tight constraints on the number of light neutrino 
species with standard-model couplings [179] through Z production in e + e _ collisions: 

IVeff = 2.984 ± 0.008. (4.1) 

1 Updated beams are presented in [83]. We have checked that using these beams (rather than the ones on 
which this analysis is based) results in less than a 0.2cr shift in parameters, as described in appendix C. 

2 The effective number of relativistic species due to three neutrinos is slightly larger than three even in the 
standard scenario due to heating caused by the injection of the entropy from the e + /e _ annihilation [e.g., 
49-52, 65, 80, 118, 148]. 
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Figure 4 . Constraints on the ACDM cosmological model using the likelihood of marginalized CMB 
bandpowers using ACT+ ACTDefl data, compared to the constraints from ACT+ACTDefl+WMAP7, 
where ACTDefl is the reconstructed deflection power spectrum from ACT data, and WMAP7 alone. 
The top six panels are fitted parameters, while the bottom two rows are derived parameters. For the 
ACT-only parameters, two cases are shown: one where we fix the optical depth t, and another where 
the scalar spectral index n s is fixed. Fixing n s tightens the constraint on the amplitude of scalar 
perturbations (the ACT+ ACTDefl, n s fixed curve lies on the ACT+ACTDefl+WMAP7 curve), but 
does not cause significant shifts in the other primary parameters, relative to leaving n s free. The 
independent CMB data sets, ACT and WMAP7, are consistent with the ACDM model. 


Relativistic species (whether neutrinos or other early relativistic species) change the 
expansion rate of the universe through their energy density and impact the perturbations in 
the early universe, affecting the damping tail of the primary CMB spectrum [9, 25, 94]. In 
the case of neutrinos, the energy density p v is lower than that of photons by a factor 

Pv/ P~i ~ (7/8)(4/H) 4 ^ 3 A r e ff , (4.2) 

where lV e ff is 3.046 in the standard ACDM model. 

Extra relativistic energy density damps the small-scale CMB power — see [94] for a 
concise recent review, and the discussions in [9, 79, 96, 97, 116, 177]. 

Figure 5 and table 2 illustrate the constraints on the !V e g- from ACT in combination 
with various probes. Previous analyses [55, 101] suggested a slight excess in the N e g. This 
preference for more damping from extra relativistic degrees of freedom is no longer present 
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Figure 5. Marginalized one-dimensional likelihood of the N e g for ACT data in combination with 
other probes. Left panel: the improvement in constraints of the ACT in combination with WMAP7 
data relative to WMAP7 data on its own. .Right panel: the N e g likelihood for WMAP7+ACT data 
and a variety of other probes. The constraints are consistent with N e g = 3.046, the value in the model 
with three neutrino species. The slight shifts in the central value can be understood in the context of 
the changes in the value of the Hubble parameter, as illustrated in figure 6. 


when analyzing the ACT 3-year data in combination with WMAP7 data. The change is 
consistent with the improved statistics of the ACT 3-year data. We find 

N eS = 2.79 ± 0.56 (WMAP7 + ACT). (4.3) 

The improvement in this value relative to the WMAP-only constraints is shown in the top 
panel of figure 6. 

The result in this analysis was obtained by imposing the consistency relation between the 
primordial helium fraction at Big Bang Nucleosynthesis (BBN) and the number of effective 
relativistic degrees of freedom [104, 167, 172, 182, see section 4.9]: 

Y P = 0.2485 + 0.0016 [(273.90 fc h 2 - 6) + 100(5 - 1)] ; 

5 - Vl + (7/43) (iVeff-3). (4.4) 

Hence, in the present analysis, the helium fraction is a determined parameter given N e g and 
the baryon density, rather than remaining fixed at the standard value of Yp — 0.24. The value 
presented in [55] was higher at N e g — 5.3 ± 1.3, as we did not impose this relation; imposing 
this constraint on the previous ACT-S data would yield a modified value of N e g — 4.3 ± 1.3, 
consistent at 1 a with the value of 3.046 expected in standard ACDM. 

The value of N e g obtained when using only the ACT-S data is closer to the value 
reported in [55]. The ACT-E data prefer a lower value for N e g, leading to a combined result 
which is ~ lcr lower than presented in [55]. 

Including the recent BAO data does not change the constraints on the relativistic species, 
while adding in the SPT data shifts the mean value to slightly higher values of N e g, but still 
consistent with the WMAP7+ACT data. The highest value for N e g comes from the addition 
of the BAO and Hubble constant data, yielding 

N eS = 3.50 ± 0.42 (WMAP7 + ACT + BAO + HST) . (4.5) 
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Figure 6. Marginalized two-dimensional (68, 95%) contours in the N e g — H a plane. The two param- 
eters are highly correlated, with larger values of H 0 leading to higher values of the effective number 
of relativistic degrees of freedom. Here BAO refers to the combined 6dF and SDSS DR9 BAO mea- 
surements, while HST is the H 0 = 73.8 ± 2.4 km s -1 Mpc -1 measurement from [153]. In contrast, 
the analysis of [34] prefers H 0 = 65.5 ± 5.5 km s -1 Mpc -1 , where the 95% error bar includes both 
systematic and statistical errors. When combining both BAO and HST measurements, A r e ff is pushed 
to higher values than in either of the individual cases, as N e g tries to reconcile the differences in the 
distance ratio r s /Dy (where r s is the co-moving sound horizon at recombination and Dy is an effec- 
tive geometric distance) between the two probes. A similar trend is seen for the matter density, Q. c h 2 . 


This is due to the degeneracy between Q c h 2 (and therefore the Hubble constant in a flat 
universe) and the relativistic species, shown in figure 6. In a flat universe, higher Ho leads to 
lower f 'l c h 2 , which increases the power on medium to small scales £ ^ 200 (as the radiation 
driving of the acoustic oscillations is reduced), leading in turn to a larger value of IV e ff 
needed to damp power in the tail of the spectrum. The mild tension between the Ho inferred 
through BAO distance measurements and the Hubble constant measurements leads to a value 
of N e g which is higher than the constraint when only the Hubble constant is added to the 
CMB data. A similar trend was seen in the recently released SPT results [95]. Smaller 
values of the Hubble constant prefer lower values of N e g [e.g. , 28, 34]. The results are 
summarized in table 3. 

In addition, the correlation between the scalar spectral index and the N e g is shown 
in figure 7. Decreasing power at small scales (through increasing iV e ff) is compensated by 
increasing the scalar spectral index, which increases small-scale power. 

Finally, we also consider an additional constraint from the recent ACT measurement 
of the skewness induced by the tSZ effect [192], The tSZ skewness signal is more sensitive 
to cr§ than any other cosmological parameter (scaling approximately as erg 1 ), allowing for 
a tight constraint with few degeneracies. Using theoretical calculations similar to those 
in [192], we find that the most significant degeneracy is with 0{,/i 2 , for which the tSZ skewness 
scales approximately as (Hf,/i 2 ) 33 . Thus, we include the erg constraint from [192] in the 
following form 

a 8 (Sl b h 2 / 0.0226) °' 3 = 0.79 ± 0.03, (4.6) 

where 0.30 = 3.3/11.1, and 11.1 is the fiducial scaling of the tSZ skewness value with erg. In 
addition to its correlation with n s , the effective number of relativistic degrees of freedom is 
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Figure 7. Marginalized two-dimensional (68, 95%) contours in the N e ff — n s plane. The two parame- 
ters are strongly correlated, with larger values of n s leading to higher values of the effective number of 
relativistic degrees of freedom. This figure is based on the 16 parameter fit of the standard cosmology 
and including variations in N e g. 



strongly correlated with erg: as N e g increases, so does ag, hence this prior lowers the effective 
number of relativistic degrees of freedom to N e g — 2.55 ± 0.40. 

The fact that ACT resolves the higher order peaks of the CMB spectrum allows for 
comparison with models that allow for departure from pure free-streaming [e.g., 40]. The 
effect on the small-scale power of a model with dark photons which are initially coupled to 
dark matter and hence only start free-streaming after they decouple implies that the phase 
shift and amplitude suppression associated with the free-streaming of radiation will not be 
uniform across all multipoles. We leave the testing of such models to future work. 


4.3 Massive neutrinos 

In the previous subsection, we estimated the number of effective relativistic species, N c g ■ If 
N e g is assumed to arise solely from neutrinos impacting the CMB, these are also assumed 
to be massless in our fiducial model. However, the CMB is sensitive to the sum of neutrino 
masses, which is related to the energy density of massive neutrinos via 


O v h 2 


Ylrrit, 

93.14 eV' 


(4.7) 


In general there will be a degeneracy between the mass of a neutrino species and increased 
relativistic degrees of freedom at early times. After the neutrinos become non-relativistic, 
neutrino free streaming washes out structure on small scales. In figure 8 we show the con- 
straints on the sum of the neutrino masses assuming N e g is fixed to 3.046. For the ACT in 
combination with WMAP7 data (and keeping the number of neutrinos fixed at N e g — 3.046), 
we find 

< 0.70eV(WMAP7 + ACT, 95%CL). (4.8) 

Adding in distance measurements through the BAO and Hubble constant prior (which breaks 
the degeneracy between and Hq) improves the constraint to 


< 0.39eV (WMAP7 + ACT + BAO + i7 0 , 95%CL) . (4.9) 
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Figure 8. Marginalized one-dimensional likelihood for the sum of the neutrino mass, Em„ in the 
case of the ACT+WMAP7 data, and with the addition of the BAO and H 0 constraint. In addition, 
we add a constraint on erg from ACT skewness measurements [192]. 


Imposing the constraint on as from the ACT skewness measurement given in eq. (4.6) yields 
< 0.30eV (WMAP7 + ACT + ACTSkewness + BAO + H 0 , 95%CL) . (4.10) 


Equations (4.8) and (4.9) show that adding the ACT data strengthens the constraints 
by factors of 1.5 and 1.8, respectively, compared to [109]. [95] find a 3cr preference for non- 
zero neutrino mass, = 0.32 ± 0.11 eV when combining WMAP7 and SPT CMB data 
with BAO, HST and SPT cluster constraints. This preference is driven by a combination 
of factors, including mild tension between the SPT and WMAP7 values of n s , between SPT 
and BAO constraints, and between the CMB and SPT cluster constraints on <7g, in the 
massless neutrino case, which reduces the error by a factor of 1.6 without changing the 
central value. The lack of evidence for non-zero neutrino mass in our analysis reflects the 
agreement between the ACT CMB, WMAP7, BAO and ACT cluster measurements when 
neutrino mass is fixed to zero. A more quantitative comparison is deferred to further work in 
light of recent Planck results. The constraint on ag from ACT clusters alone in combination 
with WMAP7 data (and BAO+Hq data) for the ACDM+Ern y model [82] is < 0.29 eV, 
for which erg = 0.802 ± 0.031. 

Figure 9 illustrates the marginalized 68% and 95% contours for the sum of the neutrino 
masses and the N e s when both are varied simultaneously. Adding in the Hubble constant 
and BAO data pushes N e g to slightly higher values, while the prior on erg from ACT tSZ 
skewness measurements lowers the effective number of degrees of freedom. In all cases the 
sum of the neutrino masses is consistent with zero. 


4.4 Early dark energy 

The small-scale damping seen in the ACT data can also be interpreted as arising from a non- 
negligible amount of dark energy at decoupling. The early dark energy (EDE) component 
may be specified through its density parameter (relative to the energy required for a flat 
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Parameter" 

ACDM 

ACDM 

ACDM 

ACDM 

ACDM 




± Nts 

+m v 

+ D e 

± a/ao 

Primary 

ioon b h 2 

2.251 ±0.047 

2.237 ±0.054 

2.242 ± 0.045 

2.257 ±0.046 

2.260 ±0.045 

ACDM 

Q c h 2 

0.114 ±0.005 

0.109 ±0.011 

0.118 ±0.006 

0.117 ±0.006 

0.115 ±0.005 


1006U 

1.040 ±0.002 

1.041 ±0.003 

1.040 ±0.002 

1.059 ±0.002 

1.060 ±0.002 


n s 

0.972 ±0.012 

0.965 ±0.020 

0.970 ±0.012 

0.977 ±0.012 

0.970 ±0.013 


T 

0.091 ±0.015 

0.091 ±0.015 

0.915 ±0.015 

0.090 ±0.015 

0.888 ±0.015 


log(10 10 A|j) 

3.19 ±0.04 

3.20 ±0.05 

3.20 ±0.05 

3.18 ±0.05 

3.20 ±0.04 

Extended 

A'eff 


2.79 ±0.56 





(eV) 



< 0.70 




fie, Wo 




< 0.03 



wo 




< -0.45 



a/ao 





1.004 ±0.005 

Derived 


0.83 ±0.03 

0.81 ±0.05 

0.77 ±0.05 

0.72 ±0.07 

0.84 ±0.03 


fi\ 

0.72 ±0.03 

0.72 ±0.03 

0.69 ± 0.04 

0.62 ±0.07 

0.73 ±0.03 



0.28 ±0.03 

0.28 ±0.03 

0.32 ±0.05 

0.38 ±0.07 

0.27 ±0.03 


Hq (km s _1 Mpc -1 ) 

70.0 ±2.4 

68.8 ±3.5 

67.2 ±3.1 

61.5 ±5.0 

71.2 ±3.6 

Secondary 

“tsz 

3.4 ±1.4 

3.4 ±1.4 

3.4 ± 1.3 

3.5 ±1.4 

3.5 ±1.4 


UkSZ 

< 8.6 

< 8.4 

< 8.5 

< 8.4 

< 8.1 


dp 

7.0 ±0.5 

7.0 ±0.5 

7.0 ±0.5 

7.2 ±0.5 

7.2 ±0.5 


d c 

5.0 ±1.0 

5.0 ±1.0 

5.0 ±1.0 

5.2 ±1.0 

5.1 ± 1.0 


d s 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 


Pc 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 


a ge 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 


a gs 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

Calibration 

y is 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.02 ±0.01 

1.01 ±0.01 


y2s 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.04 ±0.02 

1.04 ±0.02 


yle 

1.00 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.02 ±0.01 

1.02 ±0.01 


y2e 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.02 

1.00 ±0.02 


-21n£’ b 

8147 

8147 

8148 

8147 

8146 

a For one- 

tailed distributions, 

the upper 95% CL is given. 

For two-tailed 

distributions 

the 68% CL 


are shown. 

b The —2 ln«S? is based on the contribution from ACT (710 data points) and the WMAP7 likelihood (whose 
best-fit gives — 2 In =7478), including priors on the calibration parameters. The total degrees of freedom 
are therefore 8169 or 8168 if we consider the 19 or 20 parameter models. 

Table 2. ACDM and extended ‘damping tail’ parameters confidence limits from the CMB data alone: 
ACT data combined with seven-year WMAP data. In table 2 of [54], the secondary parameters, 
including calibrations, are reported assuming a best-fitting cosmological model. Here we estimate the 
calibration factors simultaneously with the ACDM parameters. 


universe) ^de( a ) and an equation of state w(a), given by [53, 191] 


^de(®) 

w(a) 


flde(O) (l g iUl °) / _ -3w 0 \ . 

O./'n'i-uO fr\\„3wn TSiefi ), 


^de(O) + fl m (0)a 3w ° 

1 din Ode (a) 

3(1 — Od e (a)) dlna 3 (a + a eq ) 


^eq 


(4.11) 


The amounts of dark energy and matter today are given by Od e (0) = OA,O m (0) = 
O c + O 5 respectively, and a eq is the scale factor at matter-radiation equality. O e is the 
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Figure 9. Marginalized two-dimensional (68, 95%) contours in the N e g — plane. The number 
of effective relativistic degrees of freedom is not strongly correlated with neutrino mass using CMB 
data alone; however, adding in complementary data sets helps tighten constraints. The star indicates 
the standard point (N c g, Em^) = (3.046,0.045 eV), where 0.045 eV is the minimum neutrino mass. 


Data set 

Aeff 

H 0 

Q c h 2 

WMAP 7+ ACT 

2.79 ±0.56 

68.8 ±3.5 

0.109 ±0.011 

WMAP7+ACT+SPT 

2.98 ±0.47 

71.2 ±3.4 

0.110 ±0.010 

WMAP7+ACT+BAO 

2.87 ±0.60 

67.7 ±3.1 

0.114 ±0.011 

WMAP7+ACT+HST 

3.18 ±0.46 

72.2 ±2.2 

0.113 ±0.011 

WMAP7+ACT+BAO+HST 

3.50 ±0.42 

71.2 ±2.1 

0.124 ±0.009 


Table 3. Constraints on N e g from the ACT data combined with seven-year WMAP data and 
other data. 

fraction of dark energy allowed at early times, while the present value of the equation of 
state of this early dark energy is expressed as wo- The scaling behavior of the equation of 
state, tracking the dominant component at each cosmic era, gives rise to a radiation-like 
component at high redshifts, enhancing the damping of the small-scale CMB tail. As the 
amount of EDE at early times tends to zero, this model approximates the standard ACDM 
cosmological model. We place a prior that w(a) > —1. This EDE model has been considered 
by many authors [29-31, 46, 47, 93, 149] showing however that current CMB observations 
combined with large scale structure data have no preference for a non-zero EDE density. We 
include dark energy clustering as discussed in [29], but we fix the dark energy sound speed 
and viscosity parameters to one and zero, respectively, as expected for a perfect fluid. Prior 
to this work, the most recent constraint on this model is from [149] who report an upper 
limit of fl e < 0.018 at 95% confidence from CMB only data, combining WMAP and SPT. 
In our analysis of the ACT data in combination with WMAP7 data and the ACT deflection 
measurement, we obtain the upper bound 

fl e < 0.025 (WMAP7 + ACT ± ACTDefl, 95%CL) , (4.12) 

with the bound on the equation of state found to be wq < —0.45 (see figure 10). We do 
not include the combination of early dark energy and N e g, as the combination of these two 
parameters is largely unconstrained using the current small-scale CMB data. 
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Figure 10. One-dimensional marginalized likelihood for the early dark energy density, fi e . Adding 
in deflection measurements (shown by the inner-most curve) tightens the constraints on the allowed 
energy density in early dark energy. 


4.5 Inflationary parameters 

Inflation provides a mechanism for the seeding of cosmological structure through small fluc- 
tuations in the early universe. We constrain the spectral index of the initial spectrum of 
scalar density fluctuations through the measurement of the high — £ tail of the angular power 
spectrum. The scalar spectrum of curvature perturbations is parameterized via [111] 

\ n s {ko)— 1+1 ln(k/ko)dn s /dliLk 

jrJ • (4-i3) 

We constrain the amplitude A^(fco), spectral index n s (k q) and ‘running’ of the spectrum 
dn s /dhik , defined at a pivot point k — 0.002 Mpc -1 . The amplitude of scalar perturbations 
is found to be 

log (10 10 A^) = 3.19 ± 0.04 (WMAP7 + ACT), (4.14) 

where the factor of 10 10 is a numerical factor included to ensure robustness to numerical 
precision errors while sampling. 

4.5.1 The scalar spectral index 

A generic prediction of inflationary models is a nearly scale invariant spectrum; any deviations 
provide powerful tests of inflationary models [8, 71, 84, 133, 134, 171]. Deviations from scale 
invariance have been tested in a variety of models and contexts [e.g., 26, 78, 123, 124, 141, 160, 
170, 178, 186, 187, 190], with the high-resolution of ACT expanding possibilities for direct 
measurements of the power spectrum [91]. The standard cosmological scenario in which no 
variation of the spectral index is allowed (i.e. dn s /d\nk = 0) provides an excellent fit to the 
data. We obtain a constraint on the scalar spectral index of 

n s = 0.972 ± 0.012 (WMAP7 + ACT). (4.15) 

When dn s /d In k is allowed to float, n s moves to slightly different values as shown in figure 11. 
When considering the ACT-S and ACT-E spectra separately (while still in combination with 
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ACDM 

ACDM 

ACDM 

ACDM 

ACDM 




+ dn s /d\nk 
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2.251 ± 0.047 
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2.248 zb 0.056 

2.236 zb 0.052 

2.246 zb 0.044 

2.279 zb 0.052 

ACDM 

Q c h 2 

0.114 ± 0.005 

0.115 ±0.006 

0.111 ±0.006 

0.113 ±0.005 

0.113 ±0.047 

0.113 ±0.005 


1006 a 

1.040 zb 0.002 

1.040 zb 0.002 

1.040 zb 0.002 

1.039 ±0.003 

1.040 zb 0.002 

1.040 zb 0.002 


n s 

0.972 zb 0.012 

0.981 ±0.031 

0.980 ±0.017 

0.967 ±0.014 

0.970 ±0.011 

0.973 ±0.012 


T 

0.091 ± 0.015 

0.092 ±0.017 

0.092 ±0.015 

0.089 ±0.015 

0.090 ±0.015 

0.092 ±0.015 


log(10 10 AO 

3.19 ±0.04 

3.19 zb 0.05 

3.16 ±0.07 

3.20 ±0.05 

3.20 ±0.04 

3.18 ±0.05 

Extended 

dn s /d\nk 


-0.004 zb 0.012 






r 



< 0.30 





Y p 




0.225 zb 0.034 




Dfc 





-0.031 ± 0.026 



G/i 






< 1.7 x 10~ 7 

Derived 


0.83 ±0.03 

0.83 ±0.03 

0.82 ± 0.03 

0.81 ± 0.03 

0.81 ±0.03 

0.82 ±0.03 


n A 

0.72 ±0.03 

0.71 ±0.03 

0.73 ± 0.03 

0.72 ±0.03 

0.63 ± 0.08 

0.73 ±0.03 



0.28 ±0.03 

0.29 zb 0.04 

0.27 ±0.03 

0.28 ±0.03 

0.40 ±0.10 

0.27 ±0.03 


Hq (km s -1 Mpc -1 ) 

70.0 zb 2.4 

69.6 ±2.9 

71.1 zb 3.1 

70.0 ±2.3 

59.8 ±8.0 

70.9 ±2.5 

Secondary 

dtSZ 

3.4 zb 1.4 

3.3 ±1.4 

3.4 zb 1.4 

3.5 ± 1.4 

3.4 zb 1.4 

3.4 ±1.3 


a kSZ 

< 8.6 

< 8.4 

< 8.5 

< 8.2 

< 8.1 

< 6.9 


dp 

7.0 ±0.5 

7.1 ±0.5 

7.0 ±0.5 

7.1 ±0.5 

7.0 ±0.5 

7.3 ±0.6 


dc 

5.0 ±1.0 

5.0 ±0.9 

5.0 ±1.0 

5.0 ±1.0 

4.9 ±0.9 

4.9 zb 1.0 


d s 

3.1 ±0.4 

3.0 ±0.5 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 


A 

2.2 zb 0.1 

2.2 zb 0.1 

2.2 zb 0.11 

2.2 zb 0.1 

2.2 zb 0.1 

2.2 zb 0.1 


cige 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 


a gs 

0.4 zb 0.2 

0.4 zb 0.2 

0.4 zb 0.2 

0.4 zb 0.2 

0.4 zb 0.2 

0.4 zb 0.2 

Calibration 

yls 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.02 zb 0.01 


yls 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.04 zb 0.02 


yle 

1.00 ±0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 

1.01 zb 0.01 


y2e 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.02 

1.00 ±0.01 

1.00 ±0.02 


-2 In if 1 

8147 

8147 

8148 

8147 

8148 c 

8149 


a For one-tailed distributions, the upper 95% CL is given. For two-tailed distributions the 68% CL 
are shown. 


b The —2 In «S? is based on the contribution from ACT (710 data points) and the WMAP7 likelihood (whose 
best-fit gives — 2 In =7478), including priors on the calibration parameters. The total degrees of freedom 
are therefore 8169 or 8168 if we consider the 19 or 20 parameter models. 

c In the case of the ACDM + Qk model, both ACT power spectrum and deflection measurements are used, 
which increases the — 2.5? value. 

Table 4. ACDM and extended single-model inflationary and damping parameters from CMB data 
alone: ACT data combined with seven-year WMAP data. 


WMAP7), we find 

71,, (ACT — E) = 0.981 ± 0.012; 

n s (ACT - S) = 0.966 ± 0.012. (4.16) 

The ACT-E data specifically prefer a higher value of the ACT-S data, which is consistent 
with n s — 0.962T0.013, the value in [55]. However, the cosmological models that describe the 
two statistically independent data sets are consistent as shown in appendix A. The difference 
in the marginalized n s is indicative of residual correlations between parameters. Combining 
the two data sets increases the value of n s relative to WMAP7 alone. It is similarly higher 
than n s — 0.9538 ± 0.0081 from WMAP7+SPT+BAO+HO reported in [173]. 

Including BAO and Ho data yields the constraint 

n s = 0.971 ± 0.009 (WMAP 7 + ACT + BAO + H 0 ) , (4.17) 
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Parameter 0 

ACDM 

ACDM 

ACDM 

ACDM 

ACDM 




± Aeff 

±m„ 

r 

+ ^ k 

Primary 

1000i,/i 2 

2.225 ±0.038 

2.271 ±0.046 

2.259 ±0.039 

2.256 ±0.040 

2.261 ±0.044 

ACDM 

tt c h 2 

0.116 ±0.003 

0.125 ±0.008 

0.115 ±0.002 

0.115 ±0.003 

0.114 ±0.005 


1006U 

1.040 ±0.002 

1.038 ±0.002 

1.040 ±0.002 

1.040 ±0.002 

1.040 ± 0.002 


n s 

0.971 ±0.009 

0.982 ±0.014 

0.974 ±0.010 

0.972 ±0.010 

0.974 ±0.011 


T 

0.090 ±0.014 

0.092 ±0.015 

0.091 ±0.014 

0.089 ±0.014 

0.092 ±0.015 


log(10 10 A|) 

3.20 ±0.03 

3.19 ±0.03 

3.19 ±0.03 

3.19 ±0.04 

3.19 ±0.04 

Extended 

iVeff 


3.50 ±0.42 





£m„ (eV) 



< 0.40 




r 




< 0.16 








-0.0020 ± 0.0047 

Derived 

<78 

0.83 ±0.02 

0.86 ±0.03 

0.79 ± 0.04 

0.83 ±0.02 

0.83 ± 0.03 


Pa 

0.71 ±0.01 

0.71 ±0.01 

0.71 ±0.01 

0.71 ±0.01 

0.72 ±0.01 



0.29 ±0.01 

0.29 ±0.01 

0.29 ±0.01 

0.29 ±0.01 

0.29 ±0.01 


H 0 (km s- 1 Mpc -1 ) 

69.3 ±1.0 

71.2 ±2.0 

68.8 ±1.0 

69.3 ±1.0 

69.0 ±1.1 

Secondary 

«tsz 

3.4 ± 1.4 

3.4 ±1.4 

3.4 ± 1.4 

3.4 ± 1.4 

3.4 ±1.5 


«kSZ 

< 8.3 

< 8.8 

< 8.5 

< 8.2 

< 8.5 


CLp 

7.0 ±0.5 

7.0 ±0.5 

7.0 ±0.5 

7.1 ±0.5 

7.0 ± 0.6 


a c 

5.0 ±1.0 

4.8 ±1.0 

5.0 ± 1.0 

5.0 ±1.0 

5.0 ±1.0 


a s 

3.1 ±0.4 

3.0 ±0.4 

3.1 ±0.4 

3.1 ±0.4 

3.1 ±0.4 


0c 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 

2.2 ±0.1 


a ge 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 

0.9 ±0.2 


a gs 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

Calibration 

yls 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 


y2s 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 

1.03 ±0.02 


yle 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.01 ±0.01 

1.00 ±0.01 


y2e 

1.00 ±0.02 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.01 


-2\n£‘ b 

8153 

8152 

8153 

8153 

8152 


“ For one-tailed distributions, the upper 95% CL is given. For two-tailed distributions the 68% CL 
are shown. 


b The —2 ln«Sf is based on the contribution from ACT (710 data points) and the WMAP7 likelihood (whose 
best-fit gives — 2 ln«S?=7478), including priors on the calibration parameters. The total degrees of freedom 
are therefore 8169 or 8168 if we consider the 19 or 20 parameter models. 

Table 5. Extended model parameters confidence limits from the ACT data combined with seven-year 
WMAP data, BAO and H 0 measurements. 


which rules out a scale- invariant n s — 1 spectrum at a significance of 3.2a. The ACT+WMAP7 
value of n s — 0.972 ± 0.012 is ~ 1.2a higher than the recently released WMAP9 result 
n s — 0.9579 ± 0.0081 [17, 90], however this constraint includes SPT data and the improved 
WMAP data. 

4.5.2 ACDM + Running Index 

In addition to the spectral index alone, we test for deviations from a perfect power law, in the 
form of running of the spectral index. Our cosmological models are run around a pivot point 
of fco — 0.002 Mpc -1 . Figure 11 shows the joint constraints on the spectral index and its 
running at a de-correlated pivot point fcg — 0.015 Mpc - * 1 , chosen to minimize the correlation 
between the two parameters [39]. We show the constraints for the ACT data combination 
with WMAP7, compared to the WMAP7 data alone. This relation between the indices at 
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Figure 11. Two-dimensional marginalized limits (68% and 95%) for the spectral index, n a , plotted 
at the pivot point k 0 = 0.015 Mpc -1 , and the running of the index dn s /d\nk, for ACT+WMAP7, 
compared to WMAP7 data alone. The data are consistent with the model dn s /d\nk = 0. 


Data set 

Recfast version 

n s 

r (95% CL). 

WMAP7 alone 

vl.4.2 

0.981 + 0.020 

< 0.32 


vl.5 

0.990 + 0.021 

< 0.39 

WMAP7 + ACT 

vl.5 

0.980 + 0.017 

< 0.30 


Table 6. Constraints on the spectral index and tensor-to-scalar ratio for different versions of the 
Recfast code. 

these two pivot points is given by: 

n s {h — 0.015 Mpc -1 ) at n s (&0 = 0.002 Mpc -1 ) 

+ 111(0.015/0.002)4^-. 

amk 

The ACT data are consistent with no running of the spectral index, with 

= -0.004 ± 0.012 (WMAP7 + ACT). 

a in K 

4.5.3 ACDM + Tensor modes 

The standard ACDM model assumes that the initial fluctuations are purely scalar modes; 
however, tensor fluctuations can also be seeded at early times, generating gravitational waves. 
While gravitational waves have not yet been detected, limits can be placed on the power of 
tensor fluctuations by constraining the so-called tensor-to-scalar ratio 

r = A 2 h (k 0 )/A 2 n (k 0 ), (4.20) 

where Af^ko) is the amplitude of primordial gravitational waves, measured at a reference 
scale ko — 0.002 Mpc -1 ; these tensor fluctuations contribute to the CMB in both temperature 
and polarization. Small-scale temperature measurements such as those from ACT help to 
constrain the amplitude r by reducing its degeneracy with the scalar spectral index, n s . 

This degeneracy arises as models with large values of r have more power from tensor 
modes at small £ and therefore require smaller values of the spectral index to compensate. The 


(4.18) 

(4.19) 
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Figure 12. Top panebTwo-dimensional marginalized distribution (68% and 95% CL) for the tensor- 
to-scalar ration r and the scalar spectral index n s . The brown contours show the constraints using 
only WMAP7 data; the blue contours show the constraints when adding the 3-year ACT data, and 
BAO and H 0 measurements, and find the constraint r < 0.16. Recfast vl.5 was used to com- 
pute the constraints. Bottom panel: one-dimensional marginalized contours for the tensor-to-scalar 
ratio for the ACT data in combination with WMAP7 (dashed curve) and for the combination of 
WMAP7+ACT+BAO-|-i7o (solid curve). 


analysis of the 300 square degree ACT 1-year data combined with WMAP7 data constrained 
r < 0.25 (95% CL) using only CMB data. In this 3-year analysis we find the limit on r has 
increased, with 

r < 0.30 (WMAP7 + ACT, 95% CL) . (4.21) 

We find an upper bound ss 15% higher for the tensor-to-scalar ratio than reported in [55]. 
This is driven by the n s — r degeneracy and the fact that the ACT 3-year constraint on n s 
is higher than the previously reported value. The curvature of the marginalized likelihood is 
shown in the bottom panel of figure 12 when only using CMB data. Once distance measure- 
ments are added the bound tightens. It is worth noting that the WMAP7-only constraints on 
the tensor-to-scalar ratio and scalar spectral index also differ from the published results [109], 
as can be seen by comparing figure 12 in this work to figure 13 in [114]. These arise from 
changes in the version of Recfast used in the WMAP7 version of this particular parameter 
combination. These differences are summarized in table 6. 

When adding in data from BAO and the Hubble constant value, we obtain 

r < 0.16 (WMAP7 + ACT + BAO+i7 0 , 95% CL) . (4.22) 

This value is generally consistent with the recently released WMAP9 results [17, 90], which 
find a tighter bound of r < 0.13 using WMAP9 data in combination with small-scale CMB 
spectra, BAO measurements and the Hubble constant measurement, of r < 0.13. 

4.6 Curvature 

One of the simplest extensions of the primary cosmological model is allowing the curvature 
of the universe to vary, through the relation 

fi*(*) = O fc (l + £) 2 , (4.23) 

where is the value of the curvature at z = 0. A negative sign for implies a closed 
universe, while > 0 describes an open universe. [70] describe how a measurement of 
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non-zero curvature places limits not only on the geometry of the universe but also on the 
space of allowed models of inflation; should positive curvature be detected at the level of 
flfc < — 10~ 4 , the current paradigm of eternal inflation would be falsified, as this scenario 
generically predicts positive values of D/.. A non-zero measurement of negative curvature 
(fifc > 0) allows one instead to place constraints on the different possible pre-inflationary 
histories and the probability measures within those frameworks. 

When combining the CMB data with measurements of the Hubble constant and BAO 
we constrain the curvature to be 

= -0.0020 ± 0.0047 (WMAP7 + ACT + BAO+H 0 ) , (4.24) 


a value consistent with a flat universe, and the previously reported WMAP7 result of 

= — 0.0023+oloo56 I 1 "]- 

4.7 ACDM + the fine structure constant a 

Given the damping of the power spectrum observed at large multipoles, one might ask if 
this damping can be explained in the context of models which allow for the variation of 
fundamental constants, such as the fine structure constant, a — e 2 /h. Dark energy is a 
possible source of variation of the fine structure constant, and hence limits on the temporal 
variation in a put limits on the equation of state at late times [e.g., 138]. The parameter a 
mediates the strength of the electromagnetic interaction, and therefore it affects the formation 
of CMB anisotropies by changing the scattering rate of the Thomson scattering process as 

f — x e n e caT, (4-25) 

where aT — 8n(a/m e ) 2 /3 is the Thomson cross section for h — c—1 and x e is the electron frac- 
tion, itself dependent on the scale factor of the universe, a{t) [e.g., 14, and references therein]. 

We quantify this physical mechanism in terms of the ratio a/ao, where ao — 1/(137.036) 
is the standard local value [132] and a is the value during the recombination process. A value 
a/ao > 1 w iH generate a shift in the acoustic peaks towards higher multipoles as well as less 
damping. [127] allowed for a variation in a and illustrated the degeneracy between such 
variations and, for example, A/ff , at early times. If variations in a arise because of some 
coupling of the photon to a scalar field [16], in general one might expect fluctuations in a 
to be both temporal and spatial [165]; however, for simplicity we ignore spatial variations 
here. In figure 13, we show the constraints on the variation in a for the ACT+WMAP7 data 
combination, which give a/ao — 1-004 ± 0.005, consistent with the standard picture. 

4.8 Lensing 

The lensing of the CMB photons by gravitational structure smears out the peaks and troughs 
of the CMB temperature spectrum, transferring power between modes. ACT reported the 
first detection of intrinsic CMB lensing [43] at 4cr; in many models this provides evidence for 
dark energy from the CMB alone [163], adding further support to the ACDM model. The 
South Pole Telescope (SPT) presented a detection of the lensing deflection power spectrum 
at 6.3cr [185]. The strength of the lensing can be simply parameterized through a lensing 
amplitude A k [32] 

^meas = theory , ( 4 . 26 ) 
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Figure 13. One-dimensional marginalized likelihood for a/ a 0 from the ACT data in combination 
with WMAP7. The constraints are consistent with no variation in a at recombination. 

where Cf <l> is the power spectrum of the lensing field <fi, and Al — 1 corresponds to the 
amount of lensing in the standard ACDM model. As this parameter multiplies the entire 
lensing power spectrum, it is defined over the same range in multipole space used for the 
theoretical temperature spectrum. We measure the lensing amplitude from the smearing of 
the temperature power spectrum as 

A l — 1.7 ± 0.38 (WMAP7 + ACT) . (4.27) 

When adding in the deflection spectrum measurement, this value is reduced to 

A l = 1.34 ± 0.23 (WMAP7 + ACT + ACTDefl) , (4.28) 

as is shown in figure 14. While this value is somewhat high, it is consistent with Al — 1 at 
the 1.5cr level, where Al — 1 is the amplitude of lensing in the ACDM scenario. We have 
tested that both the ACT-S and ACT-E patches yield the same value of Al at the 0.3a 
level, and that the constraints are independent of whether a linear or logarithmic prior on 
Al is used. When using the WMAP9 data, and an alternate binning for the ACT data, the 
central value of Al shifts to Al — 1.4 ± 0.34 (WMAP9+ACTrebin) and Al — 1.15 ± 0.21 
(WMAP9+ ACTrebin+Defl) , which is consistent with unity at the 1.2a level for the TT 
only and the 0.75a level including deflection. This alternative binning was designed to 
yield similar bandpower window functions to the ACT-S and ACT-E regions. The shift 
coming from the alternative binning procedure is ~ 0.2a for Al and less than 0.1a for other 
cosmological models. 

4.9 Primordial helium 

The conditions at decoupling are highly sensitive to the expansion rate of the universe at 
early times. Extra relativistic species change the expansion rate of the universe at a given 
temperature. An increased expansion rate cools the universe faster, and nuclear freeze-out 
occurs at an earlier time (recall that H — 1/(2 1) at early times), when the concentration of 
neutrons is higher, changing the primordial abundances of light elements including 4 He. We 
parameterize the primordial helium mass fraction as 

Y P = 2x, (4.29) 
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Figure 14. One-dimensional marginalized likelihood for the lensing parameter from WMAP7 and 
ACT temperature and deflection spectrum measurements. The black curves show the updated con- 
straints when applying an alternative binning to the ACT data (designed so that the two ACT regions 
have similar band powers), using the beams presented in [83] and when including the WMAP9 data 
with ACT. 



Figure 15. Left panel: marginalized one-dimensional constraints on the helium fraction Y p while 
fixing the number of effective relativistic degrees of freedom, N e s = 3.046. The dashed line indicates 
the standard value of Y p = 0.24. Right panel: marginalized two-dimensional 68% and 95% contours 
for the effective number of relativistic degrees of freedom and the primordial helium abundance. The 
line indicates the BBN consistency relation (given in eq. (4.4)), while the star labels the ACDM model 
with JV e fj = 3.046 and Y p = 0.24. 


where x — n n / ( n n +n p ) is the neutron abundance, n n is the number density of neutrons and n p 
the number density of protons [140]. The abundance of primordial helium is degenerate with 
N e s in a model with standard Big Bang Nucleosynthesis through eq. (4.4). The constraints 
from WMAP7 in combination with ACT data are: 

Y p = 0.225 ± 0.034 (WMAP7 + ACT), (4.30) 

where N e g is held fixed at the standard value of N e g — 3.046. The marginalized one- 
dimensional likelihood in this case is shown in the left-hand panel of figure 15. 

When constraining N e g in general, as in section 4.2 and tables 2 and 5, we impose the 
BBN consistency relation on the helium abundance. Here we also consider the case where 
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both parameters are allowed to vary freely. The constraints from this model are shown in 
the right-hand panel of figure 15, with the BBN consistency relation shown as a line in the 
plane. The model with N e g — 3.046 and Y p — 0.24 is shown as a star. 

4.10 Cosmic strings 

While cosmic strings are no longer a viable seed of all cosmic structure, should a sub-dominant 
connected gas of cosmic strings exist, they would add power to the small-scale tail of the CMB 
power spectrum. The amplitude of diffuse power is proportional to the string tension, G/x, and 
hence small-scale measurements provide a handle to constrain various string models. While 
direct searches for map linear string signatures in CMB maps provide interesting constraints 
on the string tension [66, 100], the simplest prescription for including the diffuse power is to 
model it as a network of string segments which produce loops that decay into radiation. 

We consider a Nambu string template as used in [13], modeled as a power law inversely 
proportional to multipole on especially small scales [89]. Various authors model the string 
network using different approaches [e.g., 13, 18, 22, 61, 135, 144, 184], where the equations 
of motion are solved using either Nambu, Unconnected String Model or Abelian-Higgs meth- 
ods. Other approaches allow for more freedom [60], including the mean velocity of strings, 
correlation length of the string network and small-scale structure of the strings as parameters 
in the model. 

Previous constraints on the string tension from the ACT 1-year data in combination 
with WMAP7 limited G/x < 1.6 x 10 -7 [55], while constraints from the WMAP7 data in 
combination with SPT [57] using a “Nambu-Goto” method bound the allowed string tension 
to be G/x < 1 x 1(T 5 * 7 (95% CL). In this analysis, we constrain 

G/x < 1.15 x 10“ 7 (CMB + Defl + BAO + H 0 , 95% CL) . (4.31) 

Early analyses of CMB data at (£ < 1000) in the context of string models found that the 
model with a Harrison-Zeld’ovich spectrum with n s — 1 and cosmic strings was consis- 
tent with the data [15, 21] (and even that the six-parameter model with both strings and 
a scale-invariant perturbation spectrum fit the data as well as the ACDM model without 
strings), while more recent studies have shown that this model no longer provides an ad- 
equate fit to small-scale CMB data [184]. We find that fixing the primordial power to be 
scale invariant and combining all CMB data (WMAP7+ACT+SPT+ACTDefl) with BAO 
and the Hq measurement yields a constraint on the string amplitude of q str — 0.005, or 
G/x = 7.1 x 10~ 8 (95% CL). The difference in effective y 2 values between the two models 
(n s — 1 T f/str versus n s free and no strings; both models have the same number of parameters) 
is — 21nJ2?ACDM + 2 In «Sf ns= i strings — 32, indicating an 5.7a preference for the ACDM model 
over a model with strings (using a Nambu string template) and a scale-invariant spectrum. 
Recent work by [117] note that the string amplitude is highly correlated with N e g, with N e g 
being pushed towards much higher values when both the string tension and the number of 
effective degrees of freedom are varied together. 

5 Secondary parameters 

The assumed priors placed on the secondary parameters are discussed in detail in [54] , where 

the values for assumed priors are presented and justified. We summarize and interpret our 

constraints here. 
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5.1 Point source power 


The power spectrum includes a contribution from point sources that lie below our detec- 
tion threshold. At 148 GHz and 218 GHz, both synchrotron emission from radio galaxies 
and infrared emission from dusty galaxies contribute to the observed power. Point sources 
detected in either band are masked down to a flux limit of 15 mJy. We impose a prior of 
a s — 2.9 ± 0.4 at £ — 3000 on the residual power from radio galaxies at 148 GHz from sources 
below this threshold [67], based on models from [183]. For comparison, [149] use a model 
from [48], which would predict a residual power below 15mJy that differs by about 5% from 
the prediction of the [183] model, and so lies within the lcr limit of the prior we impose. 

For the ACDM model, the error on the unresolved radio galaxy power (at 150 GHz) is 
dominated by the prior on the amplitude of power at i — 3000, with 

a s — 3.1 ± 0.4. (5.1) 

We constrain the clustered ( a c ) and Poisson ( a p ) power from CIB sources respectively 

to be 

a c *5.0 ±1.0; a p = 7.0 ± 0.5, (5.2) 

where again, the amplitudes are for templates normalized to 1 /rK 2 at £ — 3000 and 150 GHz. 

The frequency dependence of the CIB power is modeled as the square of a modified 
blackbody [see section 2.3 of 54], with effective dust temperature fixed to Td — 9.7 K. The 
effective emissivity index is constrained to be (3 C — 2.2 ± 0.1. The corresponding constraint 
on the CIB spectral index between 150 and 220 GHz is ad = 3.68 ± 0.10, where a<j is defined 
in flux density units as S„ oc v ad . This result is in good agreement with the earlier ACT 1- 
year constraint of ad — 3.69 ± 0.14 [55], and somewhat higher than, though not inconsistent 
with, ad — 3.56 ± 0.07 (or 3.45 ± 0.11, when a tSZ-CIB correlation is allowed), reported 
by SPT [149]. 

While /3 C cannot be directly associated with the emissivity of individual CIB sources, 
we expect the CIB anisotropy measured by ACT to improve dusty star-forming galaxy con- 
straints when used in conjunction with data at higher frequencies [as in 2], and through 
ongoing cross-correlation analyses. 

5.2 Thermal SZ 

The amplitude of the thermal SZ (tSZ) power spectrum is a strong function of the amplitude 
of matter fluctuations, and is very sensitive to erg [107], scaling like erf' 3 [161, 181]. We 
model the tSZ power using a template [10] from hydrodynamic simulations of cosmological 
volumes that includes radiative cooling, star formation, and AGN feedback. For a$ — 0.8 this 
template has a predicted amplitude of atsz = 5.6 at £ — 3000. In this analysis of the ACT 
3-year data we separate the thermal and kinetic SZ components [54] into two dimensionless 
spectra normalized to unity at £q — 3000 in eq. (3.3). Hence, the theoretical expectation for 
the tSZ amplitude at this — 3000 is 5.6. We constrain the tSZ amplitude to be 


a tS z — 3.4 ± 1.4. (5.3) 

Our current result corresponds to a measurement of erg = 0.751 q'q 3 and is consistent with 
our previously reported values from a t sz of as = 0.77 ± 0.05 in [55] (which assumed a a\ 
scaling), and from ACT measurements of the three point function, erg = 0.79 ± 0.03 [192]. 
Our tSZ constraint is also consistent with the SPT-reported value of atsz = 3.3 ±1.1 [151]. 
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In the ACT analysis, using only the 1-year data, [55] we modeled the SZ power spectrum 
using only a total (dimensionless) amplitude of the SZ signal (thermal plus kinetic), which 
was constrained relative to a theoretical spectrum with as — 0.8 as asz/asz,o- 8 =o.8 = 0.85 ± 
0.26 using both 148 GHz and 218 GHz frequency bands, and which can be interpreted as 
£{£ + l)/27rC'^ otSZ = 6.8 ± 2.9 /rK 2 at £ — 3000 using the [11] template. 

However, the connection between a tSZ amplitude and constraints on as depends on 
the modeling of the cluster physics. For a fixed value of as the predicted amplitude of the 
tSZ power spectrum at £ — 3000 varies by tens of percent depending on the tSZ model [e.g., 
11, 161, 181]. All the direct measurements of cluster physics are limited to high-mass, low- 
redshift clusters [e.g. 143], which do not constrain tSZ models, since around half of the 
contribution to tSZ amplitude at £ — 3000 comes from lower mass and higher redshift clusters 
and groups [e.g., 10, 181]. 

The uncertainty in the tSZ model leads to a systematic error in as, which is not in- 
cluded in the erg error bar. Combining the power spectrum measurements with higher order 
correlations and cluster number counts will provide tighter constraints on as and the tSZ 
models [23, 88]. 

5.3 tSZ- CIB correlations 

The tSZ-CIB correlation term is important in interpreting the power from the kSZ effect, 
as the frequency dependence of the tSZ-CIB and kSZ are similar across the ACT and SPT 
frequency channels [see e.g., figure 5 of 3]. A larger value of £ removes more power at 
£ ~ 3000, which can be compensated for by an increase in the kSZ amplitude. As shown 
in figure 3, £ is largely unconstrained by the ACT data, and the choice of the £ < 0.2 prior 
has little impact on the cosmological analysis. Adopting a wider prior of £ < 0.5, estimated 
as the upper limit from the modeling by [3], degrades the kSZ amplitude constraint from 
flksz < 8.6 to aksz < 9.4 (95% CL), shown by the unfilled, dashed contours in figure 16. The 
kSZ constraint from fixing £ = 0 is ciksz < 8.0, indicated by the dotted horizontal line in 
figure 16. The anti-correlation between the tSZ-CIB and kSZ power is more apparent in the 
joint fit with the SPT spectra. 

5.4 Kinematic SZ 

The kSZ power spectrum can be separated into two components. The first is a low-redshift 
component (which we refer to as the homogenous kSZ) that is sourced by perturbations in 
the free electron density after reionization with peculiar velocities; on large scales this is 
known as the Ostriker-Vishniac (OV) effect [119, 137, 174], 

The second is a high-redshift component that is sourced by fluctuations in the ionized 
fraction as well as the electron density during reionization, which we refer to as patchy 


kSZ [68, 106]. In general these two components have different shapes as a function of £ [e.g., 
126, 129, 194], 

The upper limit on the kSZ for ACT data in combination with WMAP7 data is 


When relaxing the prior on the tSZ-CIB correlation to 0 < £ < 0.5 and including SPT 
data, we obtain the bound 


aksz < 8.6 (WMAP7 + ACT, £ < 0.2, 95% CL) . 


(5.4) 


aksz < 6.9 (WMAP7 + ACT + SPT, f < 0.5, 95% CL) . 


(5.5) 
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Figure 16. Marginalized two-dimensional 68% and 95% contours for the amplitude of the kSZ 
power and the power in the tSZ-CIB correlation for the combinations of ACT+WMAP7 and 
ACT+WMAP7+SPT. The filled contours are for the prior £ < 0.2, while the unfilled contours show 
the degrading of constraints when weakening the prior to £ < 0.5. The dotted line shows the limit on 
Oksz for £ = 0. 


We approximate the homogeneous kSZ contribution at £ — 3000, aksz.hom = 1.5 ± 0.5 using 
the scaling relation from [162]; however, we do not include the uncertainty in our analysis. 
We remove this homogeneous kSZ contribution leaving only the contribution from patchy 
kSZ, aksz, patchy < 7.1. The amplitude of the patchy kSZ places constraints on the mean 
redshift [z) and duration of (A z ) reionization via the relation given by [12] as 


®kSZ, patchy — 2.03 




(5.6) 


The curves in figure 17 illustrate the constraints on z and A z derived from combinations 
of ACT, SPT, and WMAP data assuming the [12] scaling relation. The grey and black lines 
show the 68% and 95% confidence regions from large-scale WMAP 7-year CMB polariza- 
tion measurements [114]. The ACT data mildly prefer a higher kSZ amplitude than SPT, 
which leads to an increase in the kSZ upper limit for WMAP7+ACT+SPT (dashed line in 
figure 17) compared to the WMAP7+SPT case (dotted line). The choice of £ prior drives 
the constraint when SPT data are included; constraints on the tSZ-CIB correlation and thus 
kSZ power should be significantly improved with cross-correlations of the ACT maps with 
CIB-dominated Herschel/SPIRE maps [3]. 

Recently [194] obtained A z < 7.9 at 95% confidence using SPT data, although their 
analysis used a somewhat different treatment of the tSZ-CIB correlation, definition of A z 
and reionization model to that considered here. 


6 Conclusions 

We have presented parameter constraints from the ACT 3-year data in the ACT-S and ACT- 
E regions, over three seasons, at 148 GHz and 218 GHz. Using the multi-frequency likelihood 
presented in [54], the data provide constraints on both the primary cosmological model and 
a suite of secondary parameters. We confirm that the ACT data alone are fully consistent 
with the standard ACDM cosmological model. Our cosmological model is also consistent 
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Figure 17. Constraints in the z — A z plane from the patchy reionization signal at l = 3000. The 
solid curves are for the ACT data in combination with WMAP7 data, the light (dark) dotted green 
curves are for SPT+WMAP7, under the assumption of priors on the tSZ-CIB correlation of £ < 0.5 
(£ < 0.2). The light (dark) dashed blue lines show the combination for WMAP7 data, ACT and SPT, 
also under the priors of £ < 0.5 (£ < 0.2); the allowed regions are below the colored lines. The gray 
region is excluded by Gunn-Peterson [69] observations and the WMAP7 limits on z are shown by the 
dark gray (1 <t) and black (2 a) lines. 


with data from the South Pole Telescope [101, 151]. While this work was being completed, 
updated cosmological results were presented by the WMAP team [17, 90]. Our results, based 
on the ACT data in combination with WMAP7 data, are consistent with the WMAP9 results, 
which have been tightened through the addition of two years of additional data and improved 
analysis techniques. We discuss general consistency between data sets in appendix A. 

We measure a power spectrum from the thermal Sunyaev-Zel’dovich effect that is con- 
sistent with ACDM, and place an upper bound on the power from the kinematic SZ effect. 
This amplitude provides constraints on the time and duration of reionization; future cross- 
correlation studies will improve this bound. 

When combining the ACT data with data from the seven-year release from the WMAP 
satellite, we constrain a host of extended cosmological models. We see no evidence for mean 
spatial curvature, running of the scalar perturbation spectral index, or extra relativistic 
degrees of freedom beyond three light neutrinos. The amount of primordial 4 He is also 
constrained by the ACT data, and is consistent with the standard prediction from Big Bang 
Nucleosynthesis. 

The ACT data tighten the bound on the sum of light neutrino masses to 0.39 eV, using 
a measurement of <jg from the skewness of the ACT temperature map. We bound the allowed 
contribution of tensor modes to the temperature spectrum, and place limits on the amplitude 
of power coming from cosmic strings. Small-scale temperature data from ACT constrain the 
density parameter associated with early dark energy models to be less than 3%. 
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Figure 18 . The best-fit spectra from the ACT data. Left panel: the solid curve indicates the best-fit 
Ci from the 3-year ACT data set (top) and difference from the best-fit model (bottom), while the 
dashed line shows the curves from the best-fit parameters for the ACT 1-year data, as presented 
in [55]. Right panel : the best-fit power spectra for the WMAP7+ACT-E and WMAP7+ACT-S data 
sets. The bands around each best-fit spectrum are all models which lie within la of the best-fit. The 
spectra are consistent at the 4% level. At i = 2500, 4% corresponds to about 4/iK 2 , which is about 
25% larger than our error bar at i = 2500. 
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A ACT data consistency 

In order to check consistency of the current best-fit cosmological model with previous ACT 
constraints presented in [55], it is useful to consider the Ci plot derived from the best- 
fit model. This is shown in the left-hand panel of figure 18. The reduced error bars at 
larger multipoles yield tighter constraints on the positions of the fifth and sixth acoustic 
peaks, providing tighter constraints on the baryon density, fifth 2 , relative to the WMAP7- 
only constraints. This comparison is shown directly in figure 3. 


30 


JCAP10 (2013) 060 



0.0205 0.0215 0.0225 0.0235 0.0245 






3.05 


3.15 3.25 


3.35 


Figure 19. Marginalized 1-dimensional likelihoods for the ACDM cosmological parameters, com- 
puted for the ACT-E data in combination with WMAP7, the ACT-S+WMAP7 combination, and the 
constraints computed using only the ACT 1-year data, as presented in [55]. The dot-dashed curves 
show the parameters from the ACT 1-year analysis [54], while the solid and dashed lines are the 
parameters from the WMAP7+ACT-E and WMAP7+ACT-S constraints respectively. While there 
are small shifts between the two parameter sets, the best-fit spectra, shown by the curves in figure 18, 
agree to the 4% level, and the parameter sets are consistent when considering the full covariance 
between parameters. 


The marginalized one-dimensional likelihoods for the ACDM parameters are shown in 
figure 19. The ~ 0.7<r shift in 9a from the [55] release relative to the 3-year ACT data 
constraints presented here can be explained by the similar shift in parameters of the first 
1-year data maps used in [55] and [44] relative to the maps released with [56]. This shift 
was investigated thoroughly, and found to be caused by improved noise treatment and map- 
making pipeline. This shift in 8 a caused a negligible (~ 0.03cr) shift in the underlying 
cosmological parameter Da- 

The best-fit curves for the ACT-E and ACT-S spectra in combination with WMAP7 
data are shown in the right-hand panel of figure 18, while the marginalized likelihoods are also 
shown in figure 19. The ACT-E data on their own (in combination with WMAP7 data) prefer 
a higher value for 8 a and the scalar spectral index. Given the ~ 4% differences in the 

best-fit CV curves, one might ask how well the combined fit compares to the fits obtained when 
considering only the ACT-E and ACT-S subsets individually in combination with WMAP7 
data. Removing all information from WMAP7 other than a prior on r yields differences of 
< 8% between the best-fit spectra obtained using only the ACT-E and ACT-S data. 

The —2 In ££ value of the best fit to ACT-S+WMAP7 data is 444.6 relative to 446.6 
when combining ACT-S and ACT-E, for 480 data points and 10 fitted parameters, leading 
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Figure 20. Marginalized 1-dimensional likelihoods for the ACDM cosmological parameters, computed 
for the ACT-E data (without WMAP7), and the ACT-S data. The spectral index n s and optical 
depth r are held fixed, and all calibration parameters are set to their best-fit values from the ACT- 
E+WMAP7 and ACT-S+WMAP7 runs respectively. We find that the data sets are statistically 
consistent when considering the full covariance between the parameters. 


to 470 degrees of freedom. When considering WMAP7+ACT-E, the —2 In Jz? value is 227.7 
compared to 229.0 when using the ACT 3-year data set in combination with WMAP7. The 
relative difference in — 21nJz? between the south-only and equatorial-only parameters and the 
respective —2 In .£? value for the combined fit is approximately 2; hence one might expect a 
roughly 1 a shift in cosmological parameters when considering the two data sets independently. 

One can go further and test for consistency between individual seasons without includ- 
ing data from WMAP, which we show in figure 20. The roughly lcr shifts in the value of the 
baryon density and ratio of the acoustic horizon to the angular diameter distance at decou- 
pling can be understood from the fact that the data are taken from different patches on the 
sky with different noise properties. The two parameter sets are statistically independent, and 
so ACT’s internal consistency can be quantitatively checked by approximating the likelihood 
as Gaussian in the cosmological parameters. The covariance of the south and equatorial 
parameter differences is then simply the sum of the south and equatorial covariances, as re- 
ported by the MCMC chains. We excise £ from this test as it is clearly prior-driven. For the 
remaining 17 parameters, we find y 2 = 8.78 with a probability to exceed this value (PTE) of 
0.88. The kinetic and thermal SZ amplitude distributions are meaningfully affected by the 
constraint that they be non-negative; excising these two parameters in addition has only a 
minimal impact, with \ 2 = 8.14 for 15 degrees of freedom (PTE=0.84). We also find that 
no particular parameter combinations stand out as being inconsistent between the south and 
equatorial data. 
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We divide each parameter difference by its standard deviation and appropriately rescale 
the combined covariance matrix (which now is unity along the diagonal); this operation 
leaves y 2 unchanged. We rotate into the eigenvector space of the re-scaled covariance and 
find that the eigenmode distribution is consistent with a Gaussian distribution; in particular 
for both the 15- and 17-parameter case, the most discrepant mode falls within the 68% 
confidence interval expected from purely Gaussian statistics. We note that the posterior 
parameter distributions need not be Gaussian, but that Gaussianity does appear to be a 
good approximation, and that the cosmological parameters derived separately from the ACT 
south and equatorial regions are fully consistent with each other. 

For the extended models, we see a shift of ~ ler in the preferred value for N e g for 
ACT-E relative to ACT-S, but the lensing parameter Al shifts by no more than 0.2cr for the 
different patches. 

B Galactic foreground treatment 

As discussed in [54], we marginalize over residual Galactic dust allowing for two different 
levels in ACT-S and ACT-E and imposing Gaussian priors on the dust amplitudes a ge — 
0.8 ± 0.2; a gs — 0.4 ± 0.2, after masking bright dust sources in map space as reported in [42]. 
We tested that using a bigger or smaller mask has a negligible effect on all the parameters. 
Moreover we masked out in the equatorial strip a “seagull” -like structure [see figure 7 in 42], 
and again found that parameters do not move if we leave this region in the equatorial map. 
Finally we cross correlate ACT and IRIS maps [131], as described in [42], and subtract a dust 
template (rescaled by the cross-correlation coefficient derived through this cross-correlation) 
to our maps, before computing the power spectrum. We compared parameters obtained with 
or without the dust template subtraction and found the impact on parameters to be negligible. 

C Calibration and beam uncertainty 

The calibration factors of the 148 GHz and 218 GHz ACT-S and ACT-E spectra are allowed 
to vary in the cosmological analysis, yielding a 1 (2)% calibration error for the 148 (218) GHz 
spectra respectively. This error is similar to or smaller than the calibration error obtained 
from calibrating the 148 GHz maps by cross-correlating with WMAP7 maps [72], which gives 
an error of 2%. We apply a prior on the calibration parameters for the 148 GHz and 218 GHz 
spectra of 1.00 ± 0.02. We test whether removing the prior has an effect on the parameters. 
This shift is negligible for the cosmological models considered here. 

Similarly we test for a dependence of the cosmological parameters on the beam error as 
measured by planet observations, by doubling the assumed beam error and comparing this 
to the case where the beam error is ignored. We find that doubling the beam error shifts the 
peaks in the foreground parameters by < 0.1a in all cases, relative to the constraints when 
no beam error is applied. 

We present updated ACT beam parameters in [83]. With these changes, the difference 
in the window function rises up to a 2% difference on We from £ — 700 to £ — 1200 but 
then remains constant. This changes the overall calibration in the spectrum, as the relative 
calibrations are taken at £ — 700. The largest shift in parameters when using these updated 
beams (as opposed to the ones employed in this analysis) occurs in the model where N e g is a 
free parameter, and leads to a shift of < 0.2a in N e g. The overall y 2 of all models decreases 
by 8 when the updated beams are used. 
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Figure 21. ACDM parameter consistency between the full likelihood and the marginalized Ce band- 
power likelihood. The distributions are consistent between the two likelihood treatments. 

D Marginalized CMB likelihood 

To compare results to the full multi-frequency likelihood, we cross-check constraints on var- 
ious models using the simpler ‘CMB-only’ likelihood from [54], which first marginalizes over 
the secondary parameters to generate CMB band powers and a covariance matrix. This 
marginalized likelihood is then used to constrain cosmological models directly without re- 
quiring any further sampling of nuisance parameters. 

D.l ACDM model comparison 

In figure 12 of [54] we show that the secondary parameters estimated using this marginalized 
foreground likelihood and assuming the ACDM model are consistent with those obtained 
from the full multi- frequency case. There is a slight (lcr) shift in estimated kSZ amplitude, 
which shifts due to its correlation with the primary CMB; however, the constraints are still 
consistent with an upper limit on the kSZ amplitude. Figure 21 shows the marginalized 
one-dimensional likelihoods for parameters in the ACDM cosmological model. The shift 
between the constraints from the full likelihood and those from the ‘CMB-only’ likelihood 
are negligible. 

D.2 Extended model comparison 

In addition to the test for the ACDM model, we also test for the consistency between the 
full likelihood and the marginalized CMB-only likelihood for the extended models including 
the number of effective relativistic degrees of freedom, !V e ff, allowing for an estimated lensing 
amplitude, A i, including the running of the spectral index, dn s /d\nk, and variation in the 
fine-structure constant a/ao . The one-dimensional likelihoods are shown in figure 22, where 
the maximum likelihood values are consistent between methods; the shifts in the mean values 
are less then 0.05cr in all cases. Both likelihood codes are available on LAMBDA; these tests 
indicate that the CMB-only likelihood returns robust results compared to the full likelihood, 
and may be used in cases where only the primary cosmological parameters are of interest. 
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Figure 22. One-dimensional likelihoods for the extension parameters in various models. The two 
likelihoods yield consistent results. The solid and dashed lines are the same as in figure 21. 
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